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Mean-field theory is used to model polyelectrolyte adsorption and the possibility of overcompen- 
sation of charged surfaces. For charged surfaces that are also chemically attractive, the overcharging 
is large in high salt conditions, amounting to 20 — 40% of the bare surface charge. However, full 
charge inversion is not obtained in thermodynamical equilibrium for physical values of the param- 
eters. The overcharging increases with addition of salt, but does not have a simple scaling form 
with the bare surface charge. Our results indicate that more evolved explanation is needed in order 
to understand polyelectrolyte multilayer built-up. For strong polymer-repulsive surfaces, we derive 
simple scaling laws for the polyelectrolyte adsorption and overcharging. We show that the over- 
charging scales linearly with the bare surface charge, but its magnitude is very small in comparison 
to the surface charge. In contrast with the attractive surface, here the overcharging is found to 
decrease substantially with addition of salt. In the intermediate range of weak repulsive surfaces, 
the behavior with addition of salt crosses over from increasing overcharging (at low ionic strength) 
to decreasing one (at high ionic strength). Our results for all types of surfaces are supported by full 
numerical solutions of the mean-field equations. 

PACS numbers: 82.35.Gh, 82.35.Rs, 61.41. +e 



I. INTRODUCTION 

Aqueous solutions containing polyelectrolytes and 
small ions are abundant in biological systems, and have 
been the subject of extensive research in recent years. 
When such a solution is in contact with an oppositely 
charged surface, adsorption of the polyelectrolyte chains 
can occur. Theoretical descriptions of polyelectrolyte 
adsorption take into account the multitude of differ- 
ent interactions and length scales. Among others they 
include electrostatic interactions between the surface, 
monomers and salt ions, excluded- volume interactions 
between monomers and entropy considerations. Al- 
though a full description of polyelectrolytes is still lack- 
ing at present, several approaches exist and use different 
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Experimental studies |2J, |25|, |2fJ have shown that ad- 
sorbing polyelectrolytes (PEs) may carry a charge greater 
than that of the bare surface, so that the overall surface- 
polyelectrolyte complex has a charge opposite to that of 
the bare charged surface. This phenomenon is known 
as overcharging (or surface charge overcompensation) by 
the PE chains. When the overcharging is large enough to 
completely reverse the bare surface charge, the resulting 
charge surplus of the complex can be used to attract a 
second type of polyelectrolyte having an opposite charge 
to that of the first polyelectrolyte layer. Eventually, by 
repeating this process, a complex structure of alternat- 
ing layers of positively and negatively charged polyelec- 
trolytes can be formed. Experimentally, multilayers con- 
sisting of hundreds of such layers can be created [3 > 
leading the way to several interesting applications. 

A theoretical description of the PE overcharging was 
proposed in Refs. [3, |27J, |2£| based on a mean-field for- 
malism. The model relies on several approximations for a 
very dilute PE solution in contact with a charged surface 
and in theta solvent condition. In the high salt limit, the 
model predicts an exact charge inversion for indifferent 
(i.e. non-interacting) surfaces. In another work [jj, the 
scaling of the adsorption parameters was derived using a 
Flory-like free energy. This work used mean field theory, 
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but with a different type of boundary conditions. Using a 
strongly (non-electrostatic) repulsive surface for the PE, 
profiles of monomer concentration and electrostatic po- 
tential have been calculated from numerical solution of 
the mean-field equations. Several scaling laws have been 
proposed and the possibility of a weak overcharging in 
low salt condition has been demonstrated. In a related 
work |10| we presented simple scaling laws resulting from 
the same mean-field equations and for the same chemi- 
cally repulsive surfaces. In particular, we addressed the 
adsorption-depletion crossover as function of added salt. 
We showed that addition of salt eventually causes the 
polyelectrolyte to deplete from the charged surface, and 
pre-empts the high-salt adsorption regime described in 
Ref. 0. 

The present article can be regarded as a sequel of 
Ref. [lflj, offering a more complete treatment of the ad- 
sorption problem for several types of charged surfaces. In 
particular, it includes the effect of surfacc-PE chemical 
interactions, and the scaling of overcharging and adsorp- 
tion in presence of added salt. These chemical interac- 
tions are found to play a crucial role in the overcharging, 
showing that a necessary condition for the formation of 
multilayers is a chemical attraction of non-electrostatic 
origin (complcxation) between the two types of PE chains 
as well as between the PE and the charged surface. For 
chemically attractive surfaces our results deviate from 
those of Refs. |HE3)E3- We find in the same solvent and 
surface conditions that the overcharging does not reach 
a full 100% charge inversion of the bare surface charge. 
It rather depends on system parameters and never ex- 
ceeds 30-40% for physical range of parameters. Like in 
Ref. , we find an increase of the overcharging with salt 
but our numerical results do not agree with the previous 
prediction. For repulsive surfaces, several scaling laws are 
obtained and agree well with the numerically obtained 
profiles. However, due to the competition between the 
electrostatic and chemical surface interactions, the over- 
charging here is usually quite small, of the order of ~ 1% 
only. 

The article is organized as follows: the mean-field 
equations, used in the past in several other works for PE 
adsorption, are reviewed in Sec. [n] The following two 
sections treat two different types of surfaces. In Sec. IIIII 
results for PE adsorption and overcharging for attractive 
surfaces are presented and compared with previous mod- 
els. In Sec . 11 VI we derive the scaling of adsorbed PE layers 
in the case of a chemically repulsive surface. In Sec. El 



we present the adsorption and overcharging for the inter- 
mediate case of a weakly repulsive surface. In particular, 
we show the dependence of PE adsorption on the solu- 
tion ionic strength. A summary of the main results and 
future prospects are presented in Sec I VII 



II. THE MEAN-FIELD EQUATIONS AND 
SURFACE CHARACTERISTICS 

Consider an aqueous solution, containing a bulk con- 
centration of infinitely long polyelectrolyte (PE) chains, 
together with their counterions and a bulk concentra- 
tion of salt ions. Throughout this paper we assume for 
simplicity that the PE are positively charged and that 
the counterions and added salt are all monovalent. The 
mean-field equations describing such an ionic solution 
have been derived elsewhere Q and are briefly re- 
viewed here: 



v 2 v 



(2) 



where the polymer order parameter, <p(r), is the square 
root of c(r), the local monomer concentration, the 
bulk monomer concentration, ?/> the local electrostatic 
potential, c sa it the bulk salt concentration, e the dielec- 
tric constant of the aqueous solution, / the monomer 
charged fraction, e the electron charge, v the second 
virial (excluded volume) coefficient of the monomers, a 
the monomer size and — 1/fcsT the inverse of the ther- 
mal energy. Eq. (JTJ is the Poisson-Boltzmann equation 
where the salt ions, counterions and monomers are re- 
garded as sources of the electrostatic potential. Eq. @ is 
the mean-field (Edwards) equation for the polymer order 
parameter <fi(r), where the excluded- volume interaction 
between monomers and external electrostatic potential 
VK 1 ") are taken into account. 

For the case of an infinite planar wall at x = 0, Eqs. 
and can be transformed into two coupled ordinary 
differential equations, which depend only on the distance 
x from the surface. 



dx 2 



k 2 sinh£ + fc^ (e ? — rj 2 ) 



(3) 
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y0 = -W-^)+/C^ (4) 

where ( = et/j/ksT is the dimensionless (rescaled) 
electrostatic potential, rf = 4> 2 /4>l the dimension- 
less monomer concentration, kT 1 — (8irlBC sa it) the 
Debye-Hiickel screening length, due to added salt concen- 
tration, fc^ 1 = (47r? B </> 2 /) a similar decay length due 
to counterions, and Ib = e 2 /eksT the Bjerrum length. 
For water with dielectric constant e = 80 at room tem- 
perature, Ib is equal to about 7A. Note that the actual 
decay of the electrostatic potential is determined by a 
combination of salt, counterions, and polymer screening 
effects. 

The solution of Eqs. © and (JIJ requires four boundary 
conditions for the two profiles. The electrostatic poten- 
tial decays to zero at infinity, £(x — > oo) = 0. At the 
surface, we have chosen to work with a constant surface 
potential. Namely, a conducting surface with a potential 
ip = ijj s , or in rescaled variables, C(0) = — ]C»1 * The re- 
sults can be easily extended to the case of fixed surface 
charge density (which amounts to fixing the derivative of 
the potential d£/dx| x _ = — AttIbct). 

For the PE profile we also have two boundary condi- 
tions. At infinity, rj(x — ► oo) = 1, because 4>(oo) = <fib, 
has to match the bulk value. The special case of a zero 
bulk monomer concentration can be treated by taking 
4>b - *• and working directly with the non-rescaled con- 
centration, <f>(x). 

We model separately two types of surfaces. Although 
the surface always attracts the oppositely charged PE 
chains, it can be either chemically repulsive or attractive. 
In the case of a chemical repulsion between the PE chains 
and the surface, the amount of PE chains in direct con- 
tact with the surface is set to a small value rj(x = 0) = rj s . 
In the limit of a strong repulsive surface, r] s tends to zero 
and the boundary condition is r](x = 0) = 0. Because 
of the ever-present longer-range electrostatic attraction 
with the surface, PE chains accumulate in the surface 
vicinity, resulting in a positive slope dr)/dx\ x _ Q > 0. 

For chemical attractive surfaces we rely on a boundary 
condition often used 0, for neutral polymer chains: 
dij(0) /dx + r](0)/d= where d has units of length and is 
inversely proportional to the strength of non-electrostatic 
interactions of the PEs with the surface. The limit of 
d — > oo is the indifferent (non-interacting) surface limit, 
while the previous limit of a strong repulsive surface is 
obtained by a small and negative d. 



In the following sections we present our results, first for 
the chemically attractive surface and then for the chem- 
ically repulsive one. 

III. OVERCHARGING OF CHEMICALLY 
ATTRACTIVE SURFACES 

In this section we restrict the attention to surfaces that 
attract the polymer in a non-electrostatic (short range) 
fashion. For example, a possible realization can be a sys- 
tem containing PE chains with hydrophobic groups (like 
polystyrene sulfonate) that are attracted to a hydropho- 
bic surface, like the water-air interface. The short range 
attraction is modeled by a surface interaction in a similar 
way that is used extensively for neutral polymers via a 
boundary condition on the polymer order parameter at 
the surface x = 0: 

,(0) = -df (5) 

where the length d was introduced in Sec. II and is in- 
versely proportional to the strength of the surface inter- 
action. 

Close to a surface the polymer has a concentration pro- 
file given in terms of the distance from the surface x. 
Within mean-field theory the adsorbed polymer amount 
is defined with respect to the bulk concentration cf) 2 to 
be: 

p OO /" CO 

r = / dx (0 2 - 4>l) =4>l dx {ry 2 - l) , (6) 
Jo Jo 

We note that the definition of T manifests one of the 
deficiencies of mean-field theory where it is not possible 
to distinguish between chains absorbed on the surface 
and those accumulated in the surface vicinity. The latter 
will be washed away when the surface is placed in a clean 
aqueous solution and does not participate in the effective 
PE surface build-up. 

Another important quantity to be used throughout this 
paper is the overcharging parameter defined as the excess 
of PE charge over that of the bare surface charge per unit 
area, a: 

Aa = -\a\ + fT, (7) 

where a < is the induced surface charge density cal- 
culated from the fixed surface potential (C' = —AttIbO') 
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in units of e, the electron charge. Note that in terms 
of the relative overcharging parameter A<j/|aj, —1 < 
Act/ I a | < corresponds to undercharging, while a pos- 
itive A<r/\a\ > to overcharging. The special value 
Aa/\a\ = 1 indicates a full charge inversion. 

We solved numerically the mean-field coupled equa- 
tions, Eqs. Q-Q) with the electrostatic boundary condi- 
tion of a fixed surface potential, |Cs| (that has a one-to- 
one correspondence with an induced surface charge den- 
sity, a) and a non-electrostatic boundary condition from 
Eq. ©. The relative overcharging Aa/\a\ is plotted in 
Fig. 1 as function of the amount of salt. The solid line 
represents our numerical results, while the dashed one 
corresponds to a previous prediction The same sys- 
tem parameters are used for both: the limit of a dilute 
PE reservoir (4>b = 0), theta solvent conditions (v = 0), a 
strong ionic strength, and an indifferent surface taken in 
the limit of d — > oo (obtained already for d > 100 A). Fig- 
ure 1 shows an increasing dependence of the overcharging 
parameter on c sa it, but does not obey any simple scaling 
law. It varies on quite a large range of values from less 
than 10% (relative to |<r|) for c sa i t — 0.1 M to about 55% 
for c sa it — 5.5 M. We never observed for reasonable val- 
ues of the parameters a full charge inversion of 100% or 
more. 

Our results are contrasted with the approximated pre- 
diction of Ref. 6]. In the high salt limit the prediction 
reads: 



result is obtained. For smaller d, the chemical attraction 
causes a bigger overcharging. The main discrepancy be- 
tween our exact solution of the mean-field equations and 
Eq. ((HJ is in the actual limiting value for d — > oo and 
the lack of charge inversion from our results (solid line 
in Fig. 2). Note also that the difference between the two 
results cannot be explained by a constant multiplicative 
factor, and that Act/|ct| is smaller for / = 1 than for 
/ = 0.2 due to the larger electrostatic repulsion between 
charged monomers. 

To complete the presentation of the attractive surface 
we show in Fig. 3 the dependence of the overcharging 
on c S ait for several different values of surface potential, 
ranging from |£ s | — 1 to 0.2. In part (a) the PE ab- 
sorbed amount T (in units of A~ 2 ) is shown while in 
part (b) the dimensionless Acr/|a| is shown. The general 
trend is an increase of both quantities with c sa it due to 
the non-electrostatic attraction of the PE chain to the 
surface, on one hand, and the screening of the monomer- 
monomer repulsion, on the other. The overcharging is 
less than 100% unless the amount of salt is unrealisti- 
cally large. Another observation can be seen from the 
behavior of A<r/|cr|. We clearly see that Act does not de- 
pend linearly on a since the three different \Q S \ give three 
different curves (no data collapse). 

This very last result should be compared with the 
chemically repulsive surfaces which is discussed next and 
for which we find that Act ~ a. 



\a\ 3df\(r\ 

where we note that in Ref. |(| all lengths are rescaled 
with aj\f%. This prediction is shown by the dashed line 
in Fig. 1. One can see that the linear dependence on salt 
concentration is very weak and that for the entire range 
of salt the result is dominated by the first and constant 
term in Eq. (JHJ), giving an overcharging of 100%-110%. 
Although the general trend of an increased in Aer/|cr| 
appears in both results, there is neither agreement in the 
values of Acr/|a| nor in the quantitative dependence on 
Csait for a wide range of c sa it values. 

The results of Fig. 1 have been done in the limit of an 
indifferent surface modeled by d 3> k . A closer exam- 
ination of the overcharging (^-dependence is presented in 
Fig. 2. Both in part (a) for / = 0.2 and in part (b) for 
/ = 1, our result and the prediction of Eq. © show a 
similar descending trend with d. Already for d as small 
as lOOA the asymptotic results for the indifferent surface 



IV. STRONG CHEMICALLY REPULSIVE 
SURFACES 

The chemical repulsion between the PE chains and the 
surface causes the amount of PE chains in direct contact 
with the surface to diminish or even be zero for the strong 
repulsive case. The latter limit is incorporated into the 
mean-field equations by taking the boundary condition 
T} s = 0, used previously in Refs. 7, 10]. 

Our assumptions for treating the adsorbed layer are 
as follows, i) Inside the adsorption layer, the elec- 
trostatic interactions are assumed to be stronger than 
the excluded- volume interactions, ii) Using results from 
Ref. 0, we assume that the electrostatic potential de- 
cays mainly via the PE adsorption and not via the salt 
(for weak enough ionic strength), iii) Another assump- 
tion is that the electrostatic potential, Q = eip/kgT can 
be written in terms of a scaling function £ = |£ s | h(x/D), 
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where |£ s | is the surface potential and D is the adsorp- 
tion layer length scale (see also Ref. 10]). iv) The last 
assumption is that the electrostatic potential £ is low 
enough so that we can employ the linear Debye-Huckel 
approximation. 



A. First Integration of the Mean-Field Equations 



Below, for the strong repulsive case, we attempt to pro- 
duce approximate solutions to Eqs. and (|12[1 . and 
compare them to numerical calculations of Eqs. @ and 
(see also Ref. 0). 

Under the above assumptions, the excluded-volumc 
term in Eq. (|12|l can be neglected, and exp(£) and cosh(£) 
can be expanded out to second order in £, yielding: 



The salt dependence of the adsorption characteristic 
can be obtained from the first integration of Eqs. J2J and 
(0. Multiplying Eq. © by d(/dx and Eq. 10} by drj/dx, 
and then integrating both equations from x to infinity 
yields: 



2 dC. 



fir 
12 



drj 
dx 



v(f> 2 b (ij 2 



f 



r\ — dx 

dx 



dx 



(9) 



(10) 



Eq. (|10(l is then multiplied by 2k 2 n / f and subtracted 
from Eq. ©. Using (d( / dx)ij 2 dx = -rj 2 ( - 
2 f^°(dr]/dx)^r]dx we get: 



dC 

dx 



k 2 a 2 

6/ 



drj 
dx 



= k 2 (cosh C - 1) 



kl (e< 



1 



1 

2/ 



k2 mCv 2 ~ -^7V(t> 2 b k 2 m (rj 2 



1 



(11) 



which can be interpreted as the local pressure balance 
equation. The first term in the LHS of Eq. (|11|) is 
the electrostatic field pressure, and the second term is 
the pressure arising from chain elasticity. The first and 
second terms in the RHS are the ideal gas pressure of 
the salt and counterions. The third term is the pres- 
sure due to the interaction between the electrostatic field 
and the monomer concentration, and the last term is the 
excluded- volume driven pressure. 

For every segment where £ is a monotonic function of 
x, a change of variables from x to £ can be performed. 
Using drj/dx = drj/d( ■ d£/dx in Eq. lfTT|) yields: 



1 



3/ 



(v'(0f) ^) 2 = 2« 2 (coshC 



1 



2k 



2 <e< 



Cv 2 )~jv4k 2 m (^ 2 -l) 2 (12) 



3/ 



dfAVdC 

dcy v d ^ 

2fc„C (1 - 



V 



K 2 ( 2 + 

+ k 2 n ( 2 - 



(13) 



The first term on the RHS relates to the salt ions, the sec- 
ond to both the monomer and counterion concentrations 
and the third term to the counterion concentration. 



B. Scaling of Potential and Concentration Profiles 

Under the above assumptions, the dominant term 
in the RHS of Eq. i|13|l is the second one related to 
the monomer and counterion contributions. This term 
changes sign from negative to positive at r\ ~ 1. For 
1] — > (close to the surface) the negative sign of the RHS 
implies that the second (negative) term of the LHS is 
dominant: dry/d^ > y3/ /k 2 n a 2 . Finally, it can be shown 
self-consistently (presented below) that the first term on 
the LHS is of order |(s| 3 while the correction terms on 
the RHS are of order |Cs| 2 - Hence, close to the surface 
we neglect the first term in the LHS. 



k 2 n 2 

3/ 



= 2k 2 m \<;\{i- v 2 )- 
K + n 2 ) c 2 



(14) 



On the charged surface r\ = and ( — — \( s \. Using the 
mean- value theorem in the interval of interest < r\ < 1 
we get di]/d(\ y=ys ~ l/|Cr|- Substituting the scaling 
hypothesis d(/dx\ x=Q ~ \£„\/D in Eq. (|Tlfl evaluated at 
the surface yields an estimate for the layer thickness D: 



D 



a 



2M~ 



ICs- 



-1/2 



1 + I \c 1 ■ Csalt 



2m ICs 



(15) 



The first term of Eq. (|15fl retrieves the result given in 



Refs. 13,110, 12, 141. The other terms include corrections 
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due to the ionic strength of the solution. As salt is added, 
the monomer-monomer electrostatic repulsion wins over 
the electrostatic attraction of the PEs to the surface, re- 
sulting in an increase in the adsorption length D, and at 
very high salt in depletion 

By changing x to the dimensionless length x/D, ne- 
glecting the excluded-volume term and inserting the po- 
tential scaling hypothesis in Eq. the scaling form for 
the monomer concentration 4> 2 rj 2 can be derived. The 
scaling form of £ (see assumption iii above) dictates a 
similar scaling form: r\ = <j)(x)/4>b — (j) A M / $1 g{x / D), 
where 4> 2 M is the peak monomer concentration, and g is a 
scaling function normalized to one at the peak and sat- 
isfies g(0) = 0. 

To find (f>jrf, the peak concentration condition dij/d^ = 
can be inserted in Eq. (|13|l . causing the second term 
in the LHS to vanish. Using the scaling for £ and D of 
Eq. JT5J) yields: 



3|G 



2a 



alt 



47r/ B a 2 



ml 



(16) 



where factors depending on the value of h, the scaling 
function for the potential, and its derivative evaluated 
at the peak position are omitted for clarity. The above 
equation is in agreement with previous results cal- 
culated in the limit of no added-salt and negligible ef- 
fect of counterions. As the amount of salt increases, the 
monomer concentration characterized by 4> 2 M decreases, 
and for large enough amount of added salt, the peak 
monomer concentration decreases below its bulk value 
— a clear sign of depletion. 

The amount of adsorbed monomers in the adsorption 
layer is now calculated as function of ionic strength 



Wo A) 



2c sa lt + f4>\ 



4/ 



6 IC-1 3 



.(17) 



where the added subscript zero denotes the known no-salt 
limits, <p Mja = 3 |G| 2 /(4irl B a 2 ) and D = a/y/6f\Cs\, of 
the maximal monomer concentration <p 2 M and adsorption 
length D, respectively 0- For no added salt, the 
adsorbed amount scales like 



<f>M, Do ~ \Csf /2 r^l^a- 1 (18) 



adsorbed amount T decreases as can be seen from the 
negative correction term in Eq. I|17|) ■ 

Using the above assumptions, the scaling of the in- 
duced surface charge is: 



(AttIb) 



-i 

dx 



Id ICI 3/2 / 1/2 



4nl B D 



l B a 



(19) 



Comparing this scaling to Eq. 1)17(1 , we see that the scal- 
ing of a resembles that of the charge carried by the ad- 
sorbed amount fT. Subtracting the two equations shows 
that the overcharging Act — fT — a scales like a as well. 

The numerically calculated salt dependence of the ad- 
sorption length D is presented in Fig. 4a. The loca- 
tion of the concentration peak, taken as D, is shown 
as a function of the bulk salt concentration c sa it- The 
length D is shown to increase with the addition of salt, 
in agreement with Eq. l|15f) and with experimental re- 
sults [3, |3(]]. The salt dependence of the overall ad- 
sorbed amount T = <j) 2 (r/ 2 — l) dx is presented in 
Fig. 4b. The adsorbed amount is shown to decrease 
steadily with the addition of salt, in agreement with Eqs. 
(|17fl . The sharp drop in T is a sign of PE depletion in 
high salt conditions, and shows that the higher terms in 
the r(c sa it) expansion are negative as well. 

We close this section by mentioning that a more elabo- 
rated treatment of the overcharging for strongly repulsive 
surfaces is presented in Appendix A. We find two differ- 
ent scaling regimes. One for electrostatically dominated 
overcharging and the second for excluded volume dom- 
inated one. For the electrostatically dominated regime 
the overcharging follows the same scaling as the overall 
PE adsorption Act ~ T ~ fa. In the excluded volume 
dominated regime its scaling depends on the excluded 
volume parameter and yields Aero ~ / \a\ ^JIbcl 2 / (v 2 ^). 
In both cases the magnitude of overcharging is very small 
as compared to the bare a. This weak overcharging is 
not sufficient to explain PE multilayer formation for the 
repulsive surfaces considered in this section. Note that 
a very different situation exists when the surfaces are 
chemically attractive as discussed in Sec. Mil 



WEAK CHEMICALLY-REPULSIVE 
SURFACES 



When salt is added, the surface potential screening is 
obtained via the PEs and salt ions. In addition, the 



In Sec. IIVI we treated the strong chemically repelling 
surface, and concluded that the PE-surface electrostatic 
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attraction is the sole drive for the adsorption. In this 
section, we relax the assumption of strong repulsive sur- 
faces, while preserving the dominance of the electrostatic 
interactions. 

The chemical interactions are added into our model via 
the amount of PE in direct contact with the adsorbing 
surface <p(x — 0) = cj) s [or in the renormalized form r](x = 
0) = T] s ]. We note that the case of weakly repulsive 
surfaces the slope of the monomer order parameter near 
the surface must be positive drj/dx\ x=0 > 0, as explained 
in Sec. HH 

We begin by noting that the mean field equations © 
and (JIJ are invariant to translations in the coordinate 
x. Namely, the equations are invariant under a transla- 
tional transformation x — > x + l, as long as the boundary 
conditions are also transformed in the same manner, i.e. 
(f)(1) = 4> s and ((I) = — \Cs\- We note that this is a prop- 
erty of the planar geometry used in this case, and that 
for different geometries such as spherical or cylindrical 
this shift symmetry is no longer present. 

Using the above symmetry, the adsorption profile char- 
acterized by 4>(x), C( x )i with boundary conditions (f)(0) = 
4> s , C(0) = — |Cs|j can be thought of as part of a larger 
profile, satisfying £(— I) = — \Q\ an d (/>(— I) = 0. Such a 
profile, in turn, is exactly similar to the one discussed in 
the previous subsection, since the two systems are con- 
nected by the above mentioned translational transforma- 
tion. Therefore, finding the above and I from the 
given <f> s , \( s \ then enables the derivation of the adsorp- 
tion characteristics in the same manner as discussed in 
Sec. El 

This strategy is demonstrated in Fig. 5. In Fig. 5a the 
monomer concentration profiles are presented as a func- 
tion of the distance from the surface x for several surface 
parameters < (f> s < 16.2</> 6 and < |Cs| < 1.0. The 
profiles were obtained by solving Eqs. 10 and Q nu- 
merically using the boundary conditions r](0) = tf> s /(f)t>, 
C(0) = — |C| close to the surface and ((x — > oo) = 0, 
77(2; — > 00) = 1 at infinity. In Fig. 5b the same profiles 
are shown, with a translation in the surface position. As 
can be seen in Fig. 5b, all three profiles collapse on ex- 
actly the same profile, showing that all three profiles have 
the same concentration and length scales despite the dif- 
ference in \( s \ and r] s . These scales can be calculated by 
using Eqs. lfT3|l . JT^Jl, as discussed in the previous sub- 
section. 

Using the scaling relations Eqs. I|15|l and the 
above symmetry yields the following connections between 



\Cs I, I and the surface boundary conditions |£ s | and 



: <Pm9 
\C\h 



I 

~D* 



(20) 
(21) 



where D*,cf>* M are the same as Eqs. (|15|) and <|16[) when 
we change |Cs| — * IC I- Manipulation of Eq. (|21|l then 
yields: 



IC 



Ca 



h(l/D*) M-T 1 WW^m)) 



(22) 



Eq. (1221 can be solved iteratively for |£*| potential as 
a function of the surface and bulk solution parameters. 
Insertion of the result into Eq. I|21|) then gives l/D*. 

It is important to note, that the inversion of Eq. I|21|) is 
only possible in the region where <f> and x are monotonic, 
i.e. when the surface monomer concentration is lower 
than the maximal monomer concentration on the profile 
(j) 2 s < <^ ax . This condition can also be expressed by 
the condition dr]/dx\ x=0 > 0, showing that the chemical 
interactions between the surface and the PE chains must 
still be repulsive for the shift strategy to be employed. 

We now proceed to solving Eq. 1)22(1 for low values of 
<t>s/4>M- Using the boundary conditions and Eqs. 10 EJ) 
the characteristic functions h(x'), g(x') can be expanded 
in powers of x' to yield g{x') — a x' + b x 13 + . . . and 
h(x') = — 1+aix' +bix' 2 + . . . . Inserting these expansions 
into Eq. 1221) and solving to second order in <j) s yields: 



ic.r = ic.i 



Ql <t>s 
do (j>M 



(23) 



Note that <\f M is now no longer the actual maximal 
monomer concentration, but is defined in Eq. 1|16[) as the 
characteristic for the monomer concentration, where the 
real surface electrostatic potential |£ s | is used rather than 
the phantom surface potential. 

We can now turn to calculate the overall adsorption. 
Returning to the initial assumption that the monomer 
concentration profile is a part of a larger profile starting 
at x = —I, the adsorbed amount can be taken as the 
total adsorbed amount from x — — I to infinity, minus 
the amount adsorbed from x — —I to the actual surface 
at x = 0: 
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VI. DISCUSSION AND CONCLUSIONS 



Ax i 



dx i 



dx i 



) ■ (24) 



The first integral in the RHS of Eq. (|2~1J) is the overall 
adsorption of PE to the phantom surface at a; = — /, while 
the second integral is the PE adsorbed amount between 
this phantom surface and the real surface at x = 0. Using 
the expansion of g(x) and Eq. (|2(jp. we can see that the 
second integral is of third order in (j> s /<f>M, an d can be 
neglected for low enough 4> s /((>m- Using Eqs. I|17fl and 
PU)|. Eq. J2U can be evaluated as: 



,3/2 



4 7 r\/2/ s a/ 1 /2 
3ai0 s 



2c sa it 



4/^ 



2ao0M,o 



1 



2c, 



salt 



2 



IC- 



IC. 



(25) 



As expected, the adsorbed amount increases with S . 
However, the salt dependence of the surface is very dif- 
ferent from the strongly repulsive surface case. For low 
amounts of added salt and high enough <p s , the term 
combining the small ions and cj) s is stronger than the 
salt term, and the adsorbed amount increases with salt. 
When 4> s is low, the addition of salt causes the adsorbed 
amount to decrease, similar to the strongly repulsive sur- 
face case. In both cases, when the amount of salt in- 
creases to a high enough value, higher order terms be- 
come dominant, and the adsorbed amount decreases, sim- 
ilar to the infinitely repulsive case in Fig. 4b. 

The salt dependence of the adsorbed amount, T, is 
shown for several <j) s values in Fig. 6a. The adsorbed 
amount is shown to always increase with the amount of 
monomers attached to the surface 4>s- For low <fi s values, 
the adsorbed amount decreases with salt, as expected 
from Eq. i|25f) and in agreement with Fig. 4b. For higher 
<p s values, the adsorbed amount is seen to increase for low 
bulk concentrations of salt and then decrease strongly. In 
Fig. 6b the relative overcharging is plotted for the same 
<p s values, showing that the relative overcharging is still 
a very small effect, even for weakly repulsive surfaces. 
This shows that the attractive chemical interactions be- 
tween the surface and the PE chains are indeed crucial 
for the multilayer formation, and not the electrostatic 
interactions alone. 



We have presented analytical and numerical calcula- 
tions of the mean-field equations describing adsorption 
of polyelectrolytes onto charged surfaces. Three surface 
situations are discussed: chemically attractive surfaces, 
and strong and weak chemical repulsive ones. 

The strongest adsorption phenomenon is seen from nu- 
merical solution of the mean-field equations for chemi- 
cally attractive surfaces. This manifests itself in a large 
overcharging of about 40% to 55% of the bare surface 
charge for high salt conditions. On the other hand we 
did not find a full charge inversion as was pred icted ear- 
lier in Ref. Q. This means that any model [^Eif which 
tries to explain multilayer formation at equilibrium needs 
to rely on non-electrostatic complexation between the 
cationic and anionic polymer chains, beside the electro- 
static interactions. So far no simple scaling form are ob- 
tained for the PE adsorbed amount in this case. However 
the adsorbed amount (and the overcharging) are shown 
to increase with the salt amount c sa it, and to decrease 
with /, the charge fraction on the PE chain as well as 
with d, which is inversely proportional to the chemical 
interaction of the surface. 

For the case of strong chemical repulsion between the 
surface and the PE chains, we find a difference in the ef- 
fect of salt addition on the width of the adsorbed layer D 
and on the adsorbed PE amount T (Fig. 4). The width 
of the PE adsorbed layer increases with addition of salt, 
while the overall adsorbed PE amount decreases with ad- 
dition of salt. This difference results from a strong de- 
crease in the monomer concentration close to the surface 
upon addition of salt. This difference between T and D 
is in agreement with experimental results of Shubin et al 
[3o| . where silica surfaces were used to adsorb cationic 
polyacrylamide (CPAM). When salt is further added to 
the solution, the PEs stop overcompensating the charged 
surface, and we are in an under-compensation regime of 
adsorption. Previous studies ^(J showed that for even 
higher amounts of salt, scaling as c sa it ~ /|0s|, the PEs 
deplete from the charged surface. For such surfaces, the 
overcharging in most cases is found to scale like the in- 
duced surface charge density a. For weakly charged PEs 
the overcharging depends on the excluded-volume pa- 
rameter and bulk monomer concentration, and scales as 
~ fa. Very weakly charged PEs are shown to adsorb 
to the charged surface, but do not overcharge it. Our 
scaling results are in agreement with numerical calcula- 
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tions of the mean-field equations. For all PE charges, the 
overcharging of the PE with respect to a repelling sur- 
face is found to be very small, of the order of 1% of the 
bare surface charge. This naturally leads to the conclu- 
sion that the overcharging relies heavily on the chemical 
interactions between the surface and the PE chains. It 
is of much smaller importance for this type of repulsive 
surfaces than for attractive ones considered above. 

For weakly repulsive surfaces, we find that the ad- 
sorbed amount increases with the addition of salt for low 
salt amounts, and decreases for high amounts of added 
salt. This is a natural interpolation between the attrac- 
tive and repulsive surface limits. Moreover, the over- 
charging of the surface charge remains low, and for high 
amounts of added salt the PE again undercompensates 
the surface charges. 

Our results can serve as a starting point for a more 
quantitative analysis of the overcharging phenomena, and 
provide for better understanding of multilayer formation. 



In this case, the overcharging can be taken as: 

/>OC 

Aa = -\a\+ff b dx(ri 2 -l)~ 
Jo 

/>oo 

ftft / dx (rf - 1) , (Al) 

J x c 

Where the first integral is the same as Eq. (JJJ , with a be- 
ing the induced surface charge on the adsorbing surface. 
The second integral in the above equation is taken from 
x = x c to infinity, and describes the PE adsorption in the 
region where the electrostatic interaction between the PE 
and the surface is repulsive. The integral from x = to 
x c balances the surface charge almost entirely for low 
amounts of added salt, using the fact that d(/dx\ x — 
(Gauss law). 

In order to find the amount of polymers in the over- 
charging layer, we begin by examining Eq. Ijlll) at the 
potential peak, x = x c , where the first term on the LHS 
vanishes. Expanding the RHS to second order in £ with- 
out neglecting the excluded-volume term yields: 
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APPENDIX A: SCALING OF THE 
OVERCHARGING LAYER IN CHEMICALLY 
REPULSIVE SURFACES 

In order to derive scaling estimates for PE adsorption 
regime, the adsorbed PE layer is divided into two sub- 
layers around the potential peak point x = x c . The 
compensation layer is defined as i < x c , and consists 
of PEs attracted to the surface mainly by electrostat- 
ics. In contrast, in the overcharging layer (x > x c ) the 
PE chains are electrostatically repelled from the surface, 
but remain in the surface vicinity solely because of their 
chain connectivity. For low enough salt concentrations, 
the amount of charge carried by the PEs in the adsorbed 
layer is much larger than that carried by the small ions. 



a 

IT 
i 

2 



v4>t (nl 



= Kc (v 2 c - i) + 



K 
V2 

h 



(A2) 



where the values at the peak are denoted by £ c = Ci x c) 
and r) c = rj(x c )- As shown in the Appendix, n c decreases 
with addition of salt. Therefore, for a large amount of 
added salt, the LHS of Eq. I|A2(1 becomes negative (this 
is true to all orders of £ c and not only to order C% as 
shown here), while the RHS is always positive, meaning 
that there is no peak in the rescaled potential £. This 
demonstrates that surface charge overcharging can only 
occur in low enough salt conditions. 

In the overcharging layer, the previous assumptions 
made in Sec. III. A about the dominance of the elec- 
trostatic interactions are not necessarily true. Conse- 
quently, the decay of the PE concentration in this re- 
gion can be governed by either one of two interactions: 
the electrostatic or the excluded-volume repulsion be- 
tween the monomers. We consider them as two lim- 
its for overcharging, i) An electrostatically dominated 
regime, v<\>\ (77^ — l) <SC /£ c , where the excluded-volume 
term is neglected in Eq. i|A2|l . ii) The excluded-volume 
dominated regime, v<f>\ (rj 2 — l) 3> f( c - Here, the elec- 
trostatic interaction between monomers in Eq. (|A2|) is 
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neglected. In both regimes the value of drj/dx at the 
peak position is estimated from the scaling of the com- 
pensation layer to be dij/dx\ x=x ~ 4>m / {4>bD). 

1. The strongly charged regime: v4>t (??c — l) /(c 

The validity criterion of this regime is similar to the 
assumptions presented in Sec. II VI A. showing that the 
overcharging layer (x > x c ) can be thought of as an ex- 
tension of the compensation layer (0 < x < x c ). Instead 
of rederiving Act from Eq. (|A2|) , we can use Eq. (|14fl and 
the expressions for D and from Eqs. i|15|) and (|16fl . 
which yield £ c ~ \C, S \ and rj 2 ~ 4 > %[/ < l > b- The resulting 
overcharging scales like the bare surface charge: 



Act ~ fT D ~ f{4> 2 M - 4>l)D ~ |ct| . 



(A3) 



Note that the effect of added salt in this regime is similar 
to the one described in Eq. (|T7|l . Namely the added salt 
lowers the surface charge overcharging, in agreement with 
experimental [3o| and numerical [lo| results. 



The intermediate charged regime: 

1) » /Co 



In this regime, the excluded- volume interactions domi- 
nate the decay of the PE concentration. Using the regime 
validity condition and neglecting the first term in the 
RHS of Eq. I|A2() yields an expression for the monomer 
concentration at x r : 



\ 3v< 



dr) 
dx 



1,2 



Kl (A4) 



Noting that the drj/dx is continuous at 1 = x c , we 
use the compensation layer scaling estimates to find 
di]/dx\ x=x ~ 4>m /{4>bD). Substituting the latter into 
Eq. I|A4(1 and expanding to first order in both the ionic 
strength 2c sa i t + f(f> 2 — {k 2 + k^/Airls and the ratio 
of the bulk and peak monomer concentrations 4> 2 / ( / ) 'm — 
l B a 2 4>l/ \Cs\ 2 yields: 



where Eq. I|15|l and l|16|l are used for D and (j> 2 ,j, respec- 
tively. 

The overcharging is calculated from Eq. J7J where the 
characteristic length scale entering the integral is the Ed- 
wards length, £ e = a/ ^/2tv4> 2 , depending only on the 
excluded- volume interactions and not on the salt. The 
overcharging Act ~ f<j% [rj 1 — l) £ e is: 



Act 



ic s i 3/2 / 3/2 r, 2 Csalt +/0| |Cs 



ml 



2irl B a 2 4>l 
3|C,| 2 



(A6) 



In the limit of no ionic strength (c sa it = and negligible 
counterion contribution), the overcharging from Eq. (|A6|) 
scales like Acto — / |ct| ^JIbo 2 / (v 2 4> 2 ). Addition of salt 
results in a decrease of overcharging, similar to what was 
shown in the previous subsection for the strongly charged 
regime. 

By comparing the corresponding expressions for /£ c 
and v(f>l(ri 2 — 1) in both regimes, we conclude that the 
boundary between the strongly charged and intermediate 
regimes occurs at / ~ 3w |£ s | /(2nlBa 2 ) (see also Ref. Q). 
This serves as a self consistency test. 

We define the adsorption excess Ar by integrating nu- 
merically the adsorbed amount from x c to 00. Note that 
in the low salt limit, /Ar ~ Act, can be thought of as the 
overcharging parameter used earlier. It is then possible 
to make a direct comparison with the two limits discussed 
above. The adsorption excess as a function of / is pre- 
sented in Fig. 7a. For small / values, Ar scales like f 1 ^ 2 , 
in agreement with Eq. (|A6|I . while for larger values of / 
the excess adsorption Ar scales like / _1 ^ 2 , in agreement 
with Eqs. H18|) and (|A3|) . These two limiting scaling be- 
haviors are shown on Fig. 7a. In Fig. 7b, Ar is presented 
as a function of | Cs | ■ The scaling Ar ~ |Cs| 3 ^ 2 also shown 
in the figure, is in agreement with Eqs. i|18fl . (|A3|) and 
(|A6() . We note, by comparing Fig. 4b to Fig. 7, that the 
overcompensating PE are of the order of less than 1% of 
the total adsorbed amount, showing that the mean-field 
overcharging of repulsive surfaces is an extremely small 
effect. 



3 icrv 



3/2 f 1/2 



2tt ^ 



2c sa it 



2M 



2-kI b cl 2 §\ 



3|Cs 



,(A5) 



3. The undercompensation threshold 

So far we presented the case where the PE layer is 
overcompensating the surface charge, and discussed it 
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in two limits of strong and intermediate charged PEs. 
In some range of system parameters the PE charges do 
not overcompensate the surface ones. The threshold for 
having this undercompensation is briefly discussed here. 

Re-examining the validity of Eq. I|A2|I . it can be seen 
that for high enough salt the third term on the right 
dominates the largest of the first two terms when: 

2 Lbo< 

up to some numerical pre-factors. Note that the same 
inequality is valid in both limits of strong and inter- 



mediate regimes, as long as we are in high salt condi- 
tions. It should be noted, that a similar scaling rule 
was found in Refs.0, 0, for the adsorption- 
depletion crossover. Numerical results show that the 
overcharging-undercompensation transitions indeed have 
the same scaling, but differ in the constant multiplying 
the scaling result. 
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FIGURE CAPTIONS 

Fig. 1 The relative overcharging, Act/ | ct \ = (fT— |ct|)/|ct|, 
is presented as a function of the amount of added 
salt, Csait- The solid line corresponds to the numer- 
ical results (Sec. IIIlJ l. while the dashed line corre- 
sponds to the predictions from Eq. (9) in Ref. Q. 
As the surface potential, |£ s |, is fixed, ct is the nu- 
merically calculated induced surface charge. The 
numerical results show a much lower overcharging 
than the predicted ones. Furthermore, the salt de- 
pendence of the overcharging is shown to be much 
stronger. The system parameters have been chosen 
to match those of Ref. p . A dilute aqueous solu- 
tion (4>b = 0), in theta solvent conditions (v = 0). 
Other parameters are e = 80, T = 300 K, |£ s | = 1.0, 
a = 5A, / = 1, d= 100A. 

Fig. 2 The relative overcharging, Act/|ct|, is plotted 
against the chemical interaction parameter d, for 
(a) / = 0.2 and (b) / = 1. The solid line corre- 
sponds to the numerically calculated relative over- 
charging, while the dashed line corresponds to the 
theoretical predictions taken from Ref. [(J. Unlike 
the previous prediction, we do not observe a full 
charge inversion. The calculations are done in high 
salt conditions c sa i t = 1 M, corresponding to a De- 
bye length of about 3 A. Other parameters used are 
the same as in Fig. 1. Although the numerical and 
predicted profiles show a similar qualitative depen- 
dence with d, they do not coincide, and converge to 
different values at d — > oo. For practical purposes, 
d = 100 A, is already a very good approximation 
for the indifferent d — > oo surface. 

Fig. 3 (a) The numerically calculated adsorbed amount 
T = J °° dec (</> 2 — 4>1) is plotted as a function of 
the amount of added salt for three values of the 
surface potential, The solid line corresponds 

to |(s| = 1.0, the dashed-dotted line to |£ a | = 0.5 
and the dashed line to \C, S \ — 0.2. The adsorbed 
amount is shown to increase both with the amount 
of added salt and with \( tS \. Other parameters used 
are as in Fig. 1. In (b) the same results are plot- 
ted for Act/|ct| as function of the amount of added 
salt. The relative overcharging Act/|ct| is seen to de- 
crease with | C|, in contrast to T where it increases 
(see part (a)). This implies that the surface charge 
ct has a stronger dependence on |£ s | than the ad- 



sorbed amount T. Therefore, the adsorbed amount 
no longer scales with ct for attractive surfaces, in 
contrast to repulsive surfaces, Eq. I|18fl . 

Fig. 4 (a) The width of the concentration profile, D, taken 
as the peak location, is presented as a function 
of, c sa it, the added-salt concentration. The dot- 
ted line corresponds to / = 0.1, the dashed line 
to f — 0.18, the dashed-dotted line to f — 0.56 
and the solid line to / = 1. Other parameters 
used are e = 80, T = 300K, |C,| = 1, a = 5 A, 
<p 2 b = 1Q- 6 A- 3 , v = 50 A 3 . The length scale of the 
adsorption D is seen to increase with the addition 
of salt. For high enough added salt concentrations, 
the concentration peak vanishes altogether, indi- 
cating that the polymer is depleted from the sur- 
face. The adsorption-depletion crossover is denoted 
by a full circle, (b) The total adsorbed amount 
T = f °° dx (cj) 2 — ~ c m D is plotted against the 
amount of added salt. The solid line corresponds 
to f = 0.03, the dashed line to f — 0.1, the dash- 
dot to f — 0.31 and the dotted line to / = 1. 
The adsorbed amount decreases slowly with salt 
for low amounts of added salt. For high concentra- 
tions of added salt the adsorbed amount decreases 
sharply to negative values, signaling an adsorption- 
depletion transition. Other parameters used are 
the same as in (a), (reproduced from Ref. (lpf). 

Fig. 5 (a) The numerically calculated monomer concen- 
tration rj 2 (x) is plotted as a function of the dis- 
tance from the charged surface x, for several val- 
ues of |d and r/s = S /</>{,. The solid line corre- 
sponds to T] s = and \( s \ = 1.0, the dashed line 
with triangle markers to r) s = 12.55 and \C sS \ = 0.6, 
and the dashed-dotted line with square markers 
to T] s — 16.2 and \( s \ — 0.34. All profiles share 
c^it = 0.01 M, / = 1, a = 5 A, v = 10 A 3 , 
4) 2 b = 10~ 6 A- 3 , e = 80 and T = 300 K. All pro- 
files are found to have the same height in the peak 
monomer concentration, despite the difference in 
the boundary conditions, (b) the same profiles as 
in part (a), after a shift in the x — position is 
used. The solid line is not shifted, the dashed line 
with triangular markers is shifted by Ax = 2.35 A 
and the dashed-dotted line with square markers is 
shifted by Ax = 4.12 A. All three profiles show a 
data collapse on the solid line profile, corresponding 

to ici = i.o. 



Fig. 6 (a) Numerically calculated adsorbed amount of 
monomers T is plotted against the added salt con- 
centration c sa it for several values of r] s — 4> s /4>b, 
for the case of weak chemically repulsive surfaces. 
The solid line corresponds to r) s — 5, the dashed 
line corresponds to rj s — 50 and the dashed-dotted 
line to r) s = 100. All profiles share \( s \ — 1.0, 
/ = 1.0, a = 5 A, v = 10A 3 , (g = 10- 8 A~ 3 , 
e = 80 and T = 300 K. The adsorbed amount is 
seen to increase with rj s for all values of c sa it. For 
low amounts of added salt, the adsorbed amount is 
seen to increase with salt, in contrast to the rj s = 
case (strong repulsive surface) in Fig. 4b, while for 
high amounts of added salt the adsorbed amount 
decreases, in agreement with Fig. 4b. (b) The rel- 
ative overcharging Aer/ \ a\ is plotted against the 
amount of added salt for the same parameters as 
in part (a). Despite the increase in the adsorbed 
amount, the relative overcharging remains a very 
small effect. At high salt concentration, the rela- 
tive overcharging becomes negative — signaling an 
under- compensation of the surface charge. 

Fig. 7 a) Numerically calculated excess adsorption Ar, 
defined as the PE adsorbed amount from the po- 
tential peak to infinity, is plotted against /. The 
squares correspond to |( s | = 0.5, and the triangles 
to |C S | = 0.6. Both profiles share <p 2 b = 10" 6 A~ 3 , 
v = 10 2 A 3 , a = 5 A, c sa it = 0.1 mM, T = 300 K and 
e = 80. The two profiles can be fitted in the low / 
region by Ar ~ f 1 / 2 (dashed line), followed by a 
high / region where Ar ~ f~ x ^ 2 (solid line). These 
scaling results are in agreement with Eqs. (| A3|> and 
(|A6f) . b) Ar is plotted against the surface poten- 
tial \( s \. The squares correspond to / = 0.1, and 
the triangles to / = 0.3. All other parameters are 
the same as in (a). The two profiles show a scaling 



Fig. 1 Shafir and Andelman: 



of Ar~|c 



,3/2 



fitted by a solid line, in agreement 



with Eqs. I|A3(I and l|A6(l . The constant prefactors 
in the fitting lines are obtained by imposing the 
condition that the fitting line passes through the 
last numerical data point in the respective regime. 




14 

Fig. 2 Shafir and Andelman: 




Fig. 3 Shafir and Andelman: 




Fig. 4 Shafir and Andelman: 




Fig. 5 Shafir and Andelman: 
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Fig. 6 Shafir and Andelman: 




Fig. 7 Shafir and Andelman: 




